Scrambling of Hartree-Fock Levels as a universal Br ownian— Motion Process 
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We study scrambling of the Hartree-Fock single-particle levels and wave functions as electrons are 
added to an almost isolated diffusive or chaotic quantum dot with electron-electron interactions. 
We use the generic framework of the induced two-body ensembles where the randomness of the 
two-body interaction matrix elements is induced by the randomness of the eigenfunctions of the 
chaotic or diffusive single-particle Hamiltonian. After an appropriate scaling of the number of 
added electrons, the scrambling of both the HF levels and wave functions is described by a universal 
function each. These functions can be derived from a parametric random matrix process of the 
Brownian-motion type. An exception to this universality occurs when an empty level just gets 
filled, in which case scrambling is delayed by one electron. An explanation of these results is given 
in terms of a generalized Koopmans' approach. 

PACS numbers: 73.23.Hk, 05.45.Mt, 73.63.Kv, 73.23.-b 



Introduction. The single-particle levels and eigen- 
functions of a chaotic or diffusive quantum dot without 
electron-electron interactions follow the predictions of 
canonical random-matrix theory (RMT) within a band 
of ^ g levels around the Fermi energy, where g is the 
dimensionless Thouless conductance [l|, Q- In an al- 
most isolated dot, the screened Coulomb interaction V 
between electrons must be taken into account. In the 
constant-interaction (CI) model, V is approximated by 
the charging energy (a constant for a fixed number of 
electrons), and a single-particle picture and RMT can 
again be used. In this paper, we go beyond the CI model 
and use the Hartree-Fock (HF) approximation. This ap- 
proximation uses a single-particle picture and yet takes 
into account both chaotic or diffusive motion and inter- 
action effects beyond the CI model. The HF equations 
depend on the one-body density and must be solved self- 
consistently. The resulting HF single-particle levels and 
eigenfunctions depend on the number of electrons on the 
dot and, therefore, change (are "scrambled") as electrons 
are added to the dot [3, 3 ■ An experimental signature of 
scrambling is the saturation of the number of correlated 
Coulomb-blockade conductance peaks as temperature T 
increases (without scrambling, that number should rise 
linearly with T) 0. 

Chaotic or diffusive motion of the electrons induces 
randomness in the two-body matrix elements of V. For 
spinless electrons, scrambling is due to the fluctuating 
part of these matrix elements (the average yields the CI 
model) . The fluctuations depend on g and on thegeome- 
try of the dot and are, therefore, non-universal Q- Here 
we evaluate two statistical measures of scrambling and 
find them to be also non-universal. This is in accord with 
expectations based on Koopmans' limit 6J , where the HF 
wave functions are assumed to remain unchanged as an 
electron is added to the dot and where the rate of spec- 
tral scrambling depends on 5 4] . In spite of this appar- 
ent non-universality, the mechanism of scrambling does 
possess universal features. Indeed, we show that there 



exists a (non-universal) scaling of the number of added 
electrons such that our statistical measures of scrambling 
become universal. Scrambling can then be described as a 
Gaussian random process (GP) of the Brownian-motion 
type with the proper symmetry T| . An exception occurs 
when an empty HF level becomes filled. We support our 
results analytically in terms of a generalized Koopmans' 
approach. 

Approach. For a fixed number of spinless electrons, the 
induced ensembles that characterize chaotic or diffusive 
systems with interactions were classified and studied in 
Ref. |3l- Here we address their properties as electrons 
are added to the dot. We use a basis of eigenstates 
of a single-particle RMT Hamiltonian /i'"^ of orthogo- 
nal symmetry with a^ and a] the associated destruction 
and creation operators and e^ the single-particle ener- 
gies. The latter obey Wigner-Dyson statistics. The anti- 
symmetrized two-body matrix elements f^.j.; of V in 
that same basis are assumed to be Gaussian-distributed 
random variables with variances a^ characterized by a 
parameter u^. In a diffusive dot u^ (x A^/g^ where A 
is the mean single-particle level spacing. As an example 
we take the third orthogonal induced ensemble defined 
in Ref. 8] for which (j'^{vfj.^j) — Av? . The Hamiltonian 
is 
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Our numerical results are for m — AO single-particle lev- 
els and 10,000 realizations. For each realization 0} of 
the ensemble, we solve the HF equations for a range of 
numbers n of electrons and obtain the single-particle HF 
levels Ea and wave functions ipa ■ 

Unfolding the HF single-particle levels. Unfolding is 
neccessary to identify unambiguously the local spectral 
fluctuations. We calculate the ensemble average of the 
HF single-particle level density (histogram in Fig. ^ 
and fit the result with a function p(e) that continues 



smoothly across the HF gap (sohd Une). We use p(e) — 
[cl?,m)yJV2,m'^ ~ An^ia + bef. For a = 0, 6 = 1 and 
c = 1 this is just Wigner's semi-circle that corresponds 
to the choice cr^(ftr|]) = 3TO/47r^ (dashed line in Fig.P). 
With A(e) = l/p(e) we define the unfolded HF single- 
particle energies as e^ = e^l /S.{ea)- 




FIG. 1: The average HF level density p(e) for n = IQ electrons 
and u — 0.1 (histogram) is compared with an empirical fit (see 
text, solid line) and Wigner's semicircle {u — 0, dashed line). 

Scrambling. We start from a dot with n — 10 electrons 
and add 6n electrons. For every HF level a, we define two 
statistical measures for scrambling. Spectral scrambling 
is measured b^ the variance of the change of the unfolded 
level energy 



cr^(a, 6n) 
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whereas wave-function scrambling is measured by the av- 
erage of the squared wave-function overlap 
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The bar denotes the ensemble average. Complete spec- 
tral scrambling is achieved when a{a,Sn) ~ 1, i.e., when 
a is comparable to the mean level spacing of the unfolded 
spectrum. 

As |u| increases, the fluctuations of the Uj^fei'^ increase, 
too, the HF levels and wave functions should, thus, vary 
more strongly with Sn, and scrambling should increase. 
And indeed, Oa " (top panels of Fig. ^ decreases mono- 
tonically and a{a, Sn) (bottom panels) increases mono- 
tonically with |u|. In Fig.|2the levels a = 3, 5, 7 (a = 25) 
are filled (empty) in the initial dot and remain so for 
all values of Sn < 10 shown whereas levels a = 14 and 
a = 18 are initially empty and become filled for Sn = A 
and Sn = 8, respectively. Here the scrambling behavior 
differs qualitatively from that elsewhere: both cr(a, Sn) 
and Oa " remain approximately unchanged when the 
empty level a becomes filled, displaying a flat section or a 
"step" . At other values of Sn, both scrambling measures 
display their usual monotonic behavior. It is obvious 
from Fig. |21 that scrambling depends on the particular 
level and on u and is, thus, not universal. 




FIG. 2: Top panels: Squared wave-function overlap o^ ' ver- 
sus number 5n of added electrons for several HF levels a. Bot- 
tom panels: Standard deviation i7{a, 5n) of the change in the 
HF energy of level a versus Sn. "Steps" are observed across 
the HF gap. The left (right) panels correspond to an interac- 
tion parameter u — 0.05 (u — 0.1). In all cases the starting 
number of electrons is n = 10. 



Universal scrambling. For fixed n, we have shown |8J 
that the HF ensemble obeys in part universal random- 
matrix statistics. We now show that upon proper rescal- 
ing of Sn, the scrambling measures also display universal 
behavior. For reasons given below, we define the scaled 
parameter 



Snes 



Au' 



A2(e„) 



Sn 
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and replot in Fig. |21the data of Fig. |21 versus Sn^g. The 
solid lines represent the same two universal functions, 
one for the wave function overlap, and the other for the 
level variance. These are defined below. In contrast to 
Fig. 121 all curves for the same scrambling measure now 
show impressive agreement with the universal one. Uni- 
versality is violated only for those values of Sn where a 
level changes from empty to filled (steps in the scram- 
bling measures for the levels a = 14, 18). 

Gaussian random processes (GP). Our scaling equa- 
tion Q and the universal correlators can be derived from 
GP. Continuous GP depend upon an external param- 
eter X and are known to describe non-interacting chaotic 
or diffusive systems which, upon rescaling x, display uni- 
versal parametric level correlations Q- For small values 
of I a; — x' I, the random Hamiltonian hjk{x) with orthog- 




FIG. 3: Same data as in Fig. |21 but as a function of SricS 
(Eq. (2J). The solid lines are predicted from a continuous 
Gaussian process with 77 = 1. To avoid overlap, the curves 
are shifted with respect to each other. The agreement is quite 
remarkable except when an empty level becomes filled (i.e., 
for a = 14,18). 



onal symmetry obeys 



hij{x)hki{x') oc {SikSji + Sii5jk)[l ~ k\x - x'l"^ + . 



(5) 



This defines the continuous GP with exponent r] ]W^ - 
The usual parametric random-matrix models discussed 
in the literature are differentiable GP for which rj = 2. 
In contradistinction, we will show that in the HF approx- 
imation, the addition of electrons corresponds to a GP 
with rj ~ 1. This exponent is characteristic of Brownian 
motion (where time plays the role of the parameter x). 

Generalized Koopmans' approach. To understand the 
universality of scrambling we use a generalized Koop- 
mans' approach. With a denoting the single-particle 
eigenstates of the ri-electron HF Hamiltonian, that 
Hamiltonian is given by 
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where the sum is over the lowest n levels. In the same 
basis, the HF Hamiltonian for (n + 1) electrons is 






f}5 



(7) 



where p^^ is the single-particle density matrix (a projec- 
tor onto the lowest {n+l) eigenstates of the HF problem 



with (n -I- 1) electrons). We expand haj 
up to first order in v. Using psp 
and Eq. ® we find 



aj around haj 
5sf}e{n + l-(3) + 0{v) 
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w^"^ where w'^.^ 
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The approximation leading to Eq. ||SJ) is tantamount to 
ignoring the change in the HF wave functions that occurs 
when an electron is added to the dot. When restricted 
to the diagonal elements of Eq. ^ , this neglect is known 
as Koopmans' approximation. Hence, our Eq. ||SJ) can be 
considered a generalized Koopmans' relation. 

Brownian Motion. To determine the statistical prop- 
erties of the matrix w, we use the last Eq. ©. For 
small \u\, the interaction matrix elements in the HF ba- 
sis were shown [3| to be uncorrelated Gaussian random 
variables with variances that are very close to their in- 
put values in the basis used in Eq. (^. This implies 
that the Wa-y^s are uncorrelated Gaussian random vari- 
ables with mean values zero and variances a'^{way) — 
2{1 + 6aj){l — i5q„+i)(1 — 6n+i'y)u'^ . The Kronecker sym- 
bols in the last two brackets are due to antisymmetriza- 
tion and imply that the matrix w carries zeros in the 
row and column labelled {n + 1). Except for this fact 
the matrices w^"^ form a Gaussian orthogonal ensemble. 
The single-particle HF wave functions (which determine 
w*^")) are not correlated with the eigenvalues Cq . There- 
fore, the matrices w*^"^ and /i*^"-' are uncorrelated. Thus, 
for small |u| and except for effects of the HF gap, the 
set of HF Hamiltonians /i("\ /i("+^), . . . with a discrete 
parameter n corresponds approximately to a discrete or- 
thogonal GP (see Ref. [i3)- We use the fact that such 
a GP can be embedded into a continuous orthogonal GP 
h{x) by the requirement that ft.*^"^ = h{xn) for a set of 
discrete points Xn with equal spacings 5x = Xn+i — Xn. 

The exponent 77 of that GP is determined by the fact 
that the matrices w'-"' and u;("-+i) are approximately 
uncorrelated. When we use Eq. (jSJ for both matri- 
ces and approximate the HF basis for (n+l) elec- 
trons by that for n electrons, this follows from the fact 
that the matrix elements z;^„+2;cn+2 and w^„+i.„„+i 
are uncorrelated in that basis. We have checked nu- 
merically that w*^"-* and z(;("+i) remain essentially un- 
correlated even when the matrix elements Van+2-an+2 
are evaluated in the HF basis of (n -I- 1) electrons. 
We conclude that our discrete GP is characterized by 

[hiXn+2) - h{Xn+l)][h{Xn+l) - h{Xn)] « ^J^^M^^W^ « 

0. This implies that the GP must have an exponent 
of ?7 — 1. Indeed, for a GP satisfying Eq. ^ we have 
[hix^+2) - /i(x„+i)][/i(x„+i) - h{xn)] cx (2"-i-l)(<5a:)". 
That correlator vanishes if and only ii rj — 1 . 

The solid lines in Fig. |21 show that our approximation 
works even though we do not deal with a true GP (w^"^ 
carries zeros in row and column (n+l)). To see why, 
we omit on the right-hand side of Eq. © in /i*^") the 

eigenvalue e^+i and in w*-"' row and column labelled n+l 
and perform the step n — > (n+l). This procedure is 
repeated with n replaced by (n + 1), etc. We obtain a 



genuine GP except for the spacing of the eigenvalues e^ . 
These obey Wigner-Dyson statistics for a < n and for 
a > (n + 1) (see Ref. 8]) while the HF gap separating 
these two sets of levels does not. Because of the gap the 
process is not a GP and we expect less mixing between 
levels below and above the gap than would otherwise 
occur. We thus expect (on top of the step seen in the 
data) a slowing-down of the decorrelation in the vicinity 
of the Fermi level. It appears that this effect is too small 
to rise above the statistical uncertainties. 

Scaling of n. Universal correlations are obtained by 
scaling the parameter x. The scaled parameter for a GP 
with exponent rj is generally given by Xgg = D^'^x where 

D = {ea{x + Sx) — ea{x)Y 1 5x^ . We identify a; with n 
and use that relation for the smallest possible change 
i5n = 1 to define (for rj — 1) the scaled parameter ricff = 



(~(n+l) 



2(")^ 



(eq —f-a) ■ Assuming that A(eQ,) does not change 
much when one electron is added and using Koopmans' 



limit, we estimate ia 



^a«+i;an+i/^(ec). Siucc o'^ {v^n+i,cn+i) = 4w^ we ob- 
tain the simple approximate scaling given in Eq. Q). 
In general, the mean level spacing in Eq. Q depends 
(weakly) on n. Although n is a discrete parameter, the 
scaled parameter Ucg can be made as close to being con- 
tinuous as we like by choosing a sufficiently small u. 



10.35 




FIG. 4: The addition of electrons to a quantum dot is analo- 
gous to Brownian motion. Shown are three HF single-particle 
levels £a versus the number n of electrons for u = 0.1 (right) 
and u = 0.01 (left). Notice the finer scale on the left. 

Discussion. The solid lines in the top (bottom) pan- 



els of F ig. 13 display the tw o universal GP curves for 
o{6x) — \{iIj{x + 6x)\'ip{x))\^ and for a{a,Sx). These are 
computed numerically from an ry = 1 GP whose param- 
eter X is scaled according to the prescription given in 
the previous paragraph and identified with Sucg. The 
wave-function correlator can be well fit by o{Sncs) ~ 
[l+6ncs/a]^^ with a = 0.23 |2|. The steps in the data for 
a = 14 (18) at i5n = 4 (8) can be understood in the frame- 
work of the generalized Koopmans approximation ^ . In 
the HF approximation of the n-electron dot, the level 
a — (n+1) is empty. In the approximation (jS]) this level 
is not coupled by the perturbation w*^"-' to any other HF 



eigenstate 7 (we recall that w. 



7,n-f 1 



is, therefore, also an eigenstate of /i'^"+^ 

'1 



(n) 

^n-l-1,7 



0) and 



'^a=n+l ^^d 1pa=n+l 



'ya=n+l- 



In particular, 

ea=n+i ~ ea=n+i and Tp'o,=^+i ~ ^a=n+i- This is why 
the measures for both spectral and wave-function scram- 
bling remain approximately unchanged between n = a—1 
and n — a. 



When plotted versus n, the single-particle HF levels 
should exhibit Brownian-motion-like behavior since the 
addition of electrons is an ry = 1 GP even for very small 
u. This is verified in Fig. ^ where we show three HF 
levels versus n for a particular realization of the induced 
ensemble. When u is reduced from u = 0.1 (right) to u = 
0.01 (left), the behavior remains non-analytic although 
on a much finer scale. 

In conclusion, we have shown that, in a chaotic or diffu- 
sive dot with interactions, the scrambling of HF levels due 
to the addition of electrons is described by a Brownian- 
motion type process. The change of the number n of 
electrons is analogous to a discrete equally-spaced time 
series. We have used a non-universal scaling of n — > UeS 
to show that the measures of scrambling of both HF lev- 
els and HF wave functions versus Ues follow closely the 
corresponding universal correlators of an 77 = 1 Gaus- 
sian random process. The only exception occurs when an 
empty level becomes filled, causing a delay in the scram- 
bling process by one electron. 
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